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In  this  work,  several  lithium  phosphate  and  thiophosphate  materials  are  modeled  to  determine  their 
optimized  lattice  structures,  their  total  energies,  and  their  electronic  structures.  Included  in  this  study 
are  materials  characterized  by  isolated  phosphate  and  thiophosphate  groups  -  U3PS4  and  Li3P04  and 
materials  characterized  by  phosphate  and  thiophosphate  dimers  -  Li4P2S6  and  Li4P206  and  U4P2S7  and 
U4P2O7.  In  addition,  the  superionic  conducting  material  Li7P3Sn  is  also  modeled  as  are  recently  discov¬ 
ered  crystalline  argyrodite  materials  U7PS6  and  Li6PS5Cl.  A  comparison  of  Li  ion  vacancy  migration  in 
L4P2S7  and  L4P2O7  shows  the  migration  energy  barriers  in  the  thiophosphate  to  be  smaller  (less  than 
one-half)  than  in  the  phosphate. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  phosphorous  oxynitride  (LiPON)  films  were  devel¬ 
oped  at  Oak  Ridge  National  Laboratory  [1-4]  as  very  promising 
solid-state  electrolytes  for  use  in  rechargeable  batteries  and  other 
applications.  The  films  have  compositions  close  to  that  of  crys¬ 
talline  Li3P04  and  ionic  conductivities  of  10-6  Son-1.  In  previous 
work  [5-9],  we  investigated  detailed  mechanisms  for  Li  ion  migra¬ 
tion  in  idealized  models  of  LiPON  based  on  Li3P04  and  related 
phosphonitrides.  In  those  materials  we  found  that  Li  ions  can 
diffuse  by  vacancy  and  interstitial  mechanisms  with  migration 
energies  of  0.3-0.7eV.  In  the  present  work,  we  report  pre¬ 
liminary  results  for  extending  this  modeling  study  to  lithium 
thiophosphate  materials  which  have  recently  received  attention 
as  promising  candidates  for  solid-state  electrolytes  [10-19]  where 
increased  ionic  conductivity  as  large  as  10-3Scm-1  has  been 
reported. 

In  Section  2,  the  computational  methods  are  briefly  described. 
In  Section  3,  results  are  presented  for  the  optimized  structures 
of  the  materials,  comparing  corresponding  phosphates  and  thio- 
phosphates  when  possible.  In  Section  4,  heats  of  formation  are 
presented  and  compared  with  available  experimental  measure¬ 
ments.  In  Section  5,  partial  densities  of  states  are  presented 
in  order  to  develop  qualitative  pictures  of  the  bonding  proper¬ 
ties  of  the  materials.  In  Section  6,  ion  migration  energies  are 
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presented  for  one  of  the  phosphate  and  thiophosphate  materi¬ 
als.  A  summary  and  some  conclusions  are  presented  in  Section 
7. 

2.  Computational  methods 

The  computational  methods  used  in  this  work  were  the  same  as 
those  used  in  our  previous  studies  of  electrolytes  related  to  LiPON. 
In  Ref.  [9],  the  choice  of  calculational  parameters  and  the  validation 
of  the  calculational  methods  are  presented.  Briefly,  we  used  density 
functional  theory  [20,21  ]  to  treat  the  electronic  states  and  the  Born- 
Oppenheimer  approximation  to  treat  the  atomic  positions  {Ra}, 
resulting  in  a  determination  of  the  “total  energy”  E({Ra})  of  the  sys¬ 
tem.  The  local  density  approximation  (LDA)  [22]  was  used  for  the 
exchange-correlation  functional.  Most  of  the  computations  were 
carried  out  using  the  PWscf  package  [23];  while  a  few  calculations 
were  performed  using  the  abinit  [24]  and  pwpaw  [25]  packages  as 
well.  Visualizations  were  constructed  using  the  OpenDX  [26]  and 
XCrySDEN  [27]  software  packages. 

Starting  from  experimental  information  for  each  material  or  a 
related  material,  restricted  optimization  of  the  total  energy  £({RQ}) 
with  respect  to  the  atomic  positions  {Ra}  and  unit  cell  parameters 
allows  us  to  determine  stable  and  meta-stable  structures  presented 
in  Section  3  and  to  estimate  the  heat  of  formation  (AH)  presented 
in  Section  4.  For  each  meta-stable  structure,  a  qualitative  picture  of 
the  valence  state  distribution  can  be  determined  from  the  averaged 
partial  densities  of  states  (Na(E) )  presented  in  Section  5.  In  addi¬ 
tion,  migration  energies  (Em)  for  Li  ion  migration  were  estimated 
using  the  “nudged  elastic  band”  (NEB)  method  [28-30],  presented 
in  Section  6. 
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Fig.  1.  Ball  and  stick  diagrams  for  unit  cells  of  Li3PS4  (left)  and  Li3P04  (right)  in  the  Pmn2\  structure. 


3.  Simulated  crystalline  structures  [31  ] 

3.1.  Li3PS4  and  Li3P04 

U3PS4  and  Li3P04  are  characterized  by  isolated  phosphate  or 
thiophosphate  groups.  Li3P04  has  been  reported  in  two  different 
orthorhombic  crystalline  forms  -  the  (3  form  [32]  has  symmetry 
Pmn2\  (#31)  and  the  7  form  [33]  has  symmetry  Pnma  (#62).  In 
our  previous  work  [6]  we  found  that  the  (3  form  is  more  stable 
by  0.03  eV  per  formula  unit,  while  more  experimental  results  are 
reported  for  the  7  form.  The  crystal  structure  of  Li3PS4  was  reported 
by  Mercier  et  al.  [34]  to  have  an  orthorhombic  symmetry  Pnma 
(#62)  which  differs  from  7-Li3P04  structure  because  of  different 
site  occupancies,  some  of  which  are  fractional.  We  determined 
optimized  total  energies  for  Li3PS4  in  several  approximations  to 
the  Mercier  structure  and  in  the  p  and  7  structures  of  Li3P04, 
and  found  the  lowest  energy  structure  to  be  that  of  p-Li3P04 
structure.  Our  results  indicate  this  structure  to  be  more  stable  than 
the  most  stable  approximation  to  the  Mercier  structure  by  0.1  eV 
per  formula  unit  and  more  stable  than  the  7-Li3P04  structure 
by  0.2  eV  per  formula  unit.  The  comparison  of  Li3PS4  and  Li3P04 
in  the  p-Li3P04  is  illustrated  in  Fig.  1  and  the  optimized  lattice 
parameters  determined  by  our  simulations  are  compared  with 
experiment  in  Table  1.  As  reported  earlier  [6]  the  calculated  lattice 
constants  are  systematically  smaller  than  experiment  as  expected 
for  calculations  of  this  type,  while  the  relative  size  parameters  of 
the  lattice,  particularly  the  fractional  coordinates  are  well  modeled 
by  the  calculations.  The  lattice  parameters  of  Li3PS4  are  found 
to  be  roughly  25%  larger  than  those  of  Li3P04  and  the  fractional 
coordinates  are  quite  similar  for  the  two  structures. 

3.2.  LiyPSs  and  Li^PS^Cl 

There  has  been  a  number  of  recent  studies  of  lithium  thio¬ 
phosphate  materials  with  additional  sulfur  and  lithium  ions  in  a 
structure  associated  with  the  mineral  argyrodite  [35-37].  The  high 


Table  1 

Lattice  parameters  for  Li3P04  and  Li3PS4  in  the  (3-Li3P04  structure,  comparing  the 
experimental  results  (“exp”)  of  Ref.  [32]  with  the  optimized  structures  (“cal”)  found 
in  this  work  (using  the  same  crystal  conventions  as  in  Ref.  [6]). 


a  (A)  b  (A) 

c(  A) 

Li3P04  (exp)  6.12 

!  5.24 

4.86 

Li3P04  (cal)  6.01 

5.15 

4.76 

Li3PS4  (cal)  7.57 

'  6.44 

6.06 

Atom 

Fractional  coordinates 

-(*.y.  z) 

Li3P04  (exp) 

Li3 P04  (cal) 

Li3PS4  (cal) 

Li  (4b) 

(0.248,  0.328,  0.986) 

(0.247,  0.328,  0.993) 

(0.243,0.315, 1.000) 

Li  (2a) 

(0.500,  0.843,  0.989) 

(0.500,  0.838,0.992) 

(0.500,  0.849,  0.986) 

P  (2a) 

(0.000,  0.824,  0.000) 

(0.000,  0.827,  0.000) 

(0.000,  0.817,  -0.001) 

O/S  (4b) 

(0.208,  0.687,  0.896) 

(0.211,0.685,0.893) 

(0.221,0.668,0.889) 

O/S  (2a) 

(0.000,  0.105,0.900) 

(0.000,  0.112,0.896) 

(0.000,  0.115,0.884) 

O/S  (2a) 

(0.500,  0.181,0.817) 

(0.500,  0.178,0.825) 

(0.500,  0.190,  0.837) 

Li  ion  conductivity  observed  [38]  in  these  materials  is  associated 
with  high-temperature  crystalline  forms  having  multiple  fractional 
occupancy  sites.  However,  these  materials  also  have  ordered  low- 
temperature  structures  which  we  have  modeled  in  this  work.  The 
low-temperature  structure  of  Li7PS6  has  been  found  [35]  to  be 
Pna2\  (#33).  The  low-temperature  structure  of  Li6PS5Cl  has  not 
been  reported,  but  it  is  reasonable  to  assume  that  it  is  similar 
to  that  reported  [36]  for  Li6PS5I  which  forms  a  monoclinic  struc¬ 
ture  with  Cc  (#9)  symmetry.  The  calculations  optimized  the  Li7PS6 
and  Li6PS5Cl  structures  and  the  results  are  illustrated  in  Fig.  2. 
For  Li7PS6  there  are  two  unbounded  S  ions  associated  with  each 
thiophosphate  group  and  for  LigPSsCl  there  is  one  unbounded  S 
ion  and  one  unbounded  Cl  ion  for  each  thiophosphate  group.  The 
optimized  lattice  parameters  for  these  materials  are  presented  in 
Table  2.  Again  we  find  the  calculated  lattice  constants  to  be  system¬ 
atically  too  small,  while  the  fractional  coordinates  are  in  reasonably 
good  agreement  with  experiment.  In  the  case  of  Li7PS6,  the  X-ray 
data  could  not  directly  determine  the  fractional  coordinates,  but 


Fig.  2.  Ball  and  stick  diagrams  for  unit  cells  of  Li7S6  in  the  Pna2i  structure  (left)  and  LieSsl  in  the  Cc  structure  (right). 
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Table  2 

(A)  Lattice  parameters  for  low-temperature  form  of  Li7PS6  in  the  Pna2\  structure 
comparing  the  experimental  results  of  Ref.  [35]  with  the  calculational  results.  (B) 
Lattice  parameters  for  low-temperature  form  of  Li6PS5Cl  in  the  Cc  structure  com¬ 
paring  the  experimental  results  of  Ref.  [36]  for  Li6PS5I  with  the  calculational  results 
for  the  Cl  material. 


A 

a  (A) 

b(A)  c(A) 

LiyPSe  (exp) 

14.08 

6.92  9.96 

LiyPSe  (cal) 

13.82 

6.90  9.62 

Atom 

Fractional  coordinates  - 

-(x.y.  z) 

LiyPSe  (exp) 

LiyPSe  (cal) 

Li  (1) 

(0.238,  0.464,  0.998) 

(0.228,  0.533,  0.948) 

Li  (2) 

(0.771,0.416,  0.047) 

(0.777,  0.401,0.048) 

Li  (3) 

(0.294,  0.768,  0.745) 

(0.314,  0.784,  0.734) 

Li  (4) 

(0.967,  0.092,  0.088) 

(0.971,0.167,  0.097) 

Li  (5) 

(0.551,0.936,  0.928) 

(0.541,0.949,  0.946) 

Li  (6) 

(0.857,  0.021,0.365) 

(0.912,0.050,  0.409) 

Li  (7) 

(0.377,  0.830,  0.270) 

(0.376,  0.848,  0.257) 

P 

(0.124,  0.753,0.250) 

(0.118,  0.775  0.245) 

S(l) 

(0.002,0.757,  0.118) 

(0.000,  0.793,0.116) 

S  (2) 

(0.251,0.251,0.624) 

(0.262,  0.270,  0.624) 

S  (3) 

(0.379,  0.504,  0.876) 

(0.391,0.510,  0.873) 

S  (4) 

(0.881,0.505,0.873) 

(0.886,  0.468,  0.866) 

S  (5) 

(0.125,0.718,0.735) 

(0.142,0.698,0.761) 

S  (6) 

B 

(0.385,  0.740,  0.506) 

(0.372,  0.762,  0.499) 

a  (A)  b  (A) 

c(A)  p 

Li6PS5I  (exp) 

12.32  7.14 

12.45  109.5° 

Li6PS5Cl  (cal) 

11.96  6.93 

11.64  110.3° 

Atom 

Fractional  coordinates  - 

(*.y.  z) 

Li6PS5I  (exp) 

Li6PS5Cl  (cal) 

Li  (1) 

(0.985,0.305,0.371) 

(1.004,  0.363,0.408) 

Li  (2) 

(0.765,0.017,  0.885) 

(0.841,0.009,  0.878) 

Li  (3) 

(0.830,  0.515,0.385) 

(0.799,  0.533,0.373) 

Li  (4) 

(0.709,  0.740,  0.066) 

(0.684,  0.749,  0.025) 

Li  (5) 

(0.980,  0.953,0.131) 

(0.982,0.913,0.169) 

Li  (6) 

(0.997,  0.415,0.679) 

(0.976,  0.459,  0.664) 

P 

(0.747,  0.743,  0.624) 

(0.745,0.764,  0.621) 

S  (1) 

(0.691,0.983,0.672) 

(0.680,  0.995,0.177) 

S  (2) 

(0.685,0.480,  0.180) 

(0.678,0.484,  0.166) 

S  (3) 

(0.921,0.248,0.174) 

(0.927,  0.228,  0.203) 

S  (4) 

(0.692,  0.734,  0.444) 

(0.701,0.782,0.435) 

S  (5) 

(0.872,  0.248,  0.490) 

(0.871,0.270,  0.500) 

I/Cl 

(-0.004,  0.747,  0.371) 

(-0.018,  0.714,  0.354) 

Table  3 

Lattice  parameters  of  Li4P2S6  (compared  with  the  experimental  measurements  of 
Ref.  [39])  and  Li4P206  in  the  P31m  structure. 


a  (A) 

C(  A) 

xP 

zs 

Li4P2S6  (exp) 

6.070 

6.577 

0.1715 

0.3237 

0.2500 

Li4P2S6  (cal) 

5.96 

6.37 

0.1758 

0.3314 

0.2483 

Li4P206  (cal) 

4.78 

5.37 

0.1951 

0.3139 

0.2615 

they  could  be  estimated  from  a  related  material  [36].  In  the  case 
of  Li6PS5Cl,  the  calculated  results  are  compared  with  the  detailed 
measurements  for  LigPSsL  which  is  expected  to  be  similar. 

3.3.  U4P2Sq  and  L14P2OQ 

The  crystal  structure  of  U4P2S6  was  described  by  Mercier  et  al. 
[39]  as  hexagonal  P63/mcm  (#193)  with  half  occupancy  of  the  P 
(4e)  sites.  Our  electronic  structure  calculations  of  the  6  possible 
configurations  of  this  unit  cell  find  the  lowest  energy  structure  to 
be  described  by  the  P31m  (#162)  structure  which  is  a  subgroup 
of  the  original  space  group.  In  this  hexagonal  group,  the  Li  ions 
are  located  at  2c  (1/3,  2/3,  0)  and  2d  (1/3,  2/3,  1/2)  sites,  the  P 
ions  are  located  at  2e  (0,  0,  zP)  sites,  and  the  S  ions  are  located 
at  6/c  (xs,  0,  z5)  sites.  The  calculated  lattice  constants  and  frac¬ 
tional  coordinate  parameters  are  listed  in  Table  3  along  with  the 
corresponding  experimental  results.  We  were  also  able  to  simu¬ 
late  a  meta-stable  phosphate  material  with  the  same  structure; 
the  optimized  lattice  parameters  for  Li4P206  are  also  listed  in 
Table  3. 

In  constrast  to  the  other  phosphates  and  thiophosphates,  an 
interesting  characteristic  of  the  optimized  L4P2S6  and  Li4P206 
structures  is  the  presence  of  a  direct  bond  between  two  P  ions. 
In  L4P2S6  the  P-P  bond  length  is  2.24  A  which  is  10%  longer  than 
the  P-S  bond  length,  while  in  Li4P206  the  P-P  bond  length  is  2.1 0  A 
which  is  40%  longer  than  the  P-0  bond  length.  This  is  illustrated  in 
Fig.  3. 

3.4.  IA4P2S7  and  Li4P20y 

While  to  the  best  of  our  knowledge,  the  crystal  structure  of 
L4P2S7  is  unknown,  L4P2O7  was  found  to  crystallize  [40]  in  the 
P I  structure  (#2).  Our  simulations  confirm  this  structure  [9]  and 
simulations  for  the  corresponding  thiophosphate  -  L4P2S7  -  show 
that  it  has  at  least  a  meta-stable  state  in  the  same  structure. 
Table  4  summarizes  the  lattice  parameters  for  these  structures 
and  Fig.  4  shows  ball  and  stick  models  of  the  unit  cell  for  these 
structures. 


Fig.  3.  Ball  and  stick  diagrams  for  unit  cells  of  Li4P2S6  (left)  and  Li4P206  (right)  in  the  P31m  structure. 
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Table  4 

Lattice  parameters  for  Li4P207  and  Li4P2S7  in  the  PI  structure,  comparing  experimental  results  (“exp”)  of  Ref.  [40]  with  optimized  structures  (“cal”)  found  in  this  work. 


a  (A) 

b  (A) 

c  (A)  a 

P 

Y 

Li4P207  (exp) 

8.56 

7.11 

5.19  111.4° 

90.0° 

103.1° 

Li4P207  (cal) 

8.40 

7.01 

5.13  111.7° 

90.0° 

103.6° 

Li4P2S7  (cal) 

10.67 

8.81 

5.80  111.1° 

89.9° 

90.8° 

Atom 

Fractional  coordinates 

-(x,y,  z) 

Li4P207(exp) 

Li4P207  (cal) 

Li4P2S7  (cal) 

Li  (1) 

(0.791,0.673,0.644) 

(0.788,0.677,  0.653) 

(0.854,  0.647,  0.476) 

Li  (2) 

(0.452,  0.755,  0.702) 

(0.444,  0.745,  0.699) 

(0.489,  0.751,0.689) 

Li  (3) 

(0.754,  0.412,0.016) 

(0.763,  0.422,  0.027) 

(0.679,  0.478,  -0.133) 

Li  (4) 

(0.988,0.854,  0.110) 

(0.987,  0.855,0.106) 

(1.012,0.880,  0.149) 

P(l) 

(0.846,  0.222,  0.438) 

(0.850,  0.227,  0.444) 

(0.801,0.216,  0.450) 

P  (2) 

(0.631,-0.214,  0.214) 

(0.626,  -0.216,  0.215) 

(0.685,  -0.189,  0.254) 

0/S(l) 

(0.684,  0.038,  0.308) 

(0.680,  0.042,0.315) 

(0.676,  0.040,  0.203) 

0/s  (2) 

(0.875,0.259,  0.741) 

(0.882,0.268,0.752) 

(0.796,  0.231,0.800) 

O/S  (3) 

(0.794,  0.398,  0.390) 

(0.798,  0.407,  0.394) 

(0.713,0.413,0.417) 

O/S  (4) 

(0.982,0.151,0.270) 

(0.985,0.149,  0.268) 

(0.975,0.173,0.301) 

O/S  (5) 

(0.622,  -0.261,0.476) 

(0.618,  -0.264,  0.482) 

(0.689,  -0.158,0.612) 

O/S  (6) 

(0.472,  -0.270,  0.047) 

(0.463,  -0.276,  0.045) 

(0.529,  -0.299,  0.071) 

O/S  (7) 

(0.763,  -0.297,  0.041) 

(0.763,  -0.296,  0.041) 

(0.839,  -0.309,  0.080) 

3.5.  Li7P3Sn 

The  crystal  structure  of  the  superionic  conductor  Li7P3Sn  was 
recently  analyzed  by  Yumane  et  al.  [13]  who  found  it  has  P\  sym¬ 
metry  (#2).  The  measured  and  calculated  lattice  parameters  are 
listed  in  Table  5.  It  is  interesting  that  the  calculated  value  of  the 
b  lattice  parameter  is  larger  than  experiment  while  the  a  and  c 
lattice  parameters  follow  the  usual  trend  of  being  smaller  than 
experiment.  The  calculated  fractional  coordinates  are  generally 
in  good  agreement  with  the  experimental  values.  Yumane  et  al. 
have  shown  that  this  structure  is  composed  of  an  ordered  arrange¬ 
ment  of  U4P2S7  and  Li3S4  groups  similar  to  the  structures  seen  in 
Figs.  4  and  1. 

4.  Heats  of  formation 

In  order  to  assess  the  chemical  stability  of  the  materials,  we 
estimated  the  heats  of  formation  relative  to  their  decomposition  to 
elemental  materials  in  their  standard  states  as  defined  in  the  CRC 
Handbook  [41],  using  the  methods  described  in  our  earlier  work 
[9].  In  order  to  extend  the  analysis  to  sulfur  containing  materi¬ 
als  we  needed  to  estimate  the  equilibrium  energy  of  the  standard 
state  of  elemental  S,  which  is  the  orthorhombic  form  (a-Sg)  [42] 
having  the  structure  Fddd  (#70).  In  order  to  extend  the  analysis 
to  Cl-containing  materials,  we  needed  to  estimate  the  energy  of 
the  standard  state  of  elemental  Cl,  which  is  molecular  Cl2.  Instead 


of  evaluating  this  energy  directly,  we  estimated  it  from  the  total 
energies  of  LiCl  in  the  rocksalt  structure,  UCIO4  in  the  Pnma  (#62) 
structure  [43],  and  PC15  in  the  P4/n  (#85)  structure  [44]  and  the 
heat  of  formation  data  values  for  these  materials  given  in  the  CRC 
Handbook  [41]. 

Some  results  are  given  in  Table  6  where  a  few  results  are 
repeated  from  Ref.  [9].  From  our  previous  experience  and  from  the 
accuracy  of  finding  the  standard  energy  of  elemental  chlorine,  we 
expect  that  the  error  of  the  calculation  to  be  ±0.5  eV  although  the 
relative  error  between  materials  of  similar  composition  is  expected 
to  be  considerably  less. 

In  addition  to  the  phosphate  and  thiophosphate  materials,  the 
reference  materials  Li20,  Li2S,  LiCl,  LiC104,  and  PCI5  are  are  also 
listed. 

From  these  results,  we  can  make  some  comments  on  relative 
stabilities  of  these  materials.  For  example,  the  calculations  suggest 
that  Li4P2S6  is  more  stable  than  L4P2S7  in  the  sense  that  in  their 
standard  states: 

AH(Li4P2S6)  +  AH(S)  =  AH(Li4P2S7)  -  0.84  eV.  (1) 

In  fact  there  are  literature  reports  of  glassy  Li4P2S7,  but  we  know 
of  no  reports  of  crystalline  Li4P2S7.  By  constrast,  for  the  analogous 
phosphate  materials,  we  find 

AH(Li4P206)  +  1aH(02)  =  AH(Li4P207)  +  4.21  eV,  (2) 


Fig.  4.  Ball  and  stick  diagrams  for  unit  cells  of  Li4P2S7  (left)  and  Li4P207  (right)  in  the  PI  structure. 
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Table  5 

Lattice  parameters  for  Li7P3Sii  in  the  PI  structure,  comparing  experimental  results 
(“exp”)  of  Ref.  [13]  with  optimized  structures  (“cal”)  found  in  this  work. 


a  (A)  b( A)  c  (A) 

a 

P  y 

Li7P3Sn 

(exp) 

12.50  6.03  12.53 

102.8° 

113.2°  74.5° 

Li7P3Sn 

(cal) 

12.00  6.16  12.24 

102.4° 

113.8°  72.1° 

Fractional  coordinates  - 

(x,y,  z) 

Atom 

Li7P3Sn  (cal) 

Li7P3Sn  (cal) 

Li  (1) 

(0.669,  0.132,  0.092) 

(0.669,  0.070,  0.094) 

Li  (2) 

(0.638,  0.402,  0.866) 

(0.600,  0.549,  0.880) 

Li  (3) 

(0.216,  0.530,  0.740) 

(0.199,  0.489,  0.707) 

Li  (4) 

(-0.058,  0.786,  0.270) 

(-0.030,  0.837,  0.251) 

Li  (5) 

(0.359,  0.591,0.337) 

(0.299,  0.771,0.298) 

Li  (6) 

(0.637,  0.734,  0.498) 

(0.680,  0.648,  0.508) 

Li  (7) 

(0.139,  0.294,  0.134) 

(0.129,  0.268,0.116) 

P(l) 

(0.794,  0.033,  0.465) 

(0.780,  0.054,  0.452) 

P  (2) 

(0.493,  0.038,  0.297) 

(0.472,-0.007,  0.273) 

P  (3) 

(0.840,  0.270,  0.059) 

(0.830,  0.318,  0.057) 

S  (1) 

(0.861,0.236,  0.620) 

(0.822,  0.289,  0.598) 

S  (2) 

(0.824,  0.105,0.334) 

(0.811,0.118,  0.314) 

S  (3) 

(0.840,  -0.317,  0.468) 

(0.853,  -0.271,0.488) 

S  (4) 

(0.618,0.160,  0.453) 

(0.583,0.170,  0.425) 

S  (5) 

(0.331,0.227,  0.288) 

(0.304,  0.152,0.294) 

S  (6) 

(0.509,  0.080,  0.151) 

(0.478,  0.047,  0.120) 

S  (7) 

(0.509,  -0.300,  0.307) 

(0.522,  -0.343,  0.302) 

S  (8) 

(0.841,0.330,  -0.097 

(0.820,  0.394,  -0.102) 

S  (9) 

(0.828,  -0.073,  0.046) 

(0.865,  -0.028,  0.057) 

S(10) 

(0.678,  0.487,  0.064) 

(0.663,0.474,  0.081) 

S  (11) 

(0.977,  0.351,0.188) 

(0.968,  0.429,  0.203) 

suggesting  that  Li4P206  is  significantly  unstable  relative  to  Li4P207. 
In  fact  we  know  of  no  literature  reports  of  observations  of  L4P2O6 
crystals. 

For  the  argyrodite  material,  the  calculations  also  suggest  that 
Li7PS6  is  unstable  relative  to  decomposition  into  Li3PS4  and  Li2S  in 
the  sense  that 

AH(Li7PS6)  =  AH(Li3PS4)  +  2 AH(Li2S)  +  0.32  eV.  (3) 

In  this  estimate,  we  used  the  (3-Li3P04  structure  to  evaluate 
AH(Li3PS4)  and  the  result  is  within  our  expected  calculation  error. 
The  material  with  chlorine  is  also  marginally  unstable  relative  to 
decomposition  into  Li3PS4,  Li2S,  and  LiCl  as  shown  in  the  equation 

AH(Li6PS5Cl)  =  AH(Li3PS4)  +  AH(Li2S)  +  AH(LiCl)  +  0.38  eV.  (4) 


Table  6 

Heats  of  formation  per  formula  unit  calculated  for  the  listed  reference  and  lithium 
phosphate  and  thiophosphate  materials.  The  structure  is  described  in  terms  of  the 
space  group  using  the  notation  of  Ref.  [31  ].  When  available,  experimental  values  of 
AH  from  Ref.  [41  ]  are  also  listed. 


Material 

Structure 

A Hcal  (eV) 

AHeXp  (eV) 

Li20 

Fm3m  (#225) 

-6.13 

-6.20 

|3-Li3P04 

Pmn2\  (#31) 

-21.31 

y-Li3P04 

Pnma  (#62) 

-21.28 

-21.72 

Li4P207 

pi (#2) 

-34.10 

Li4P206 

P31m  (#162) 

-29.89 

Li2S 

Fm3m  (#225) 

-4.26 

-4.57 

P-Li3PS4 

Pmn2\  (#31) 

-8.32 

y-Li3PS4 

Pnma  (#62) 

-8.12 

Li7PS6 

Pna2i  (#33) 

-16.51 

Li4P2S7 

PI (#2) 

-11.51 

Li4P2S6 

P31m  (#162) 

-12.35 

Li7P3Sn 

pi (#2) 

-19.89 

LiCl 

Fm3m  (#225) 

-3.80 

-4.23 

LiC104 

Pnma  (#62) 

-4.49 

-3.95 

PC15 

P4/n  (#85) 

-4.08 

-4.60 

Li6PS5Cl 

Cc  (#9) 

-16.00 

t - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 


-8  -6  -4  -2  0  2  4  6  8 

E  (eV) 


Fig.  5.  Partial  densities  of  states  for  Li3PS4  and  Li3P04. 

The  calculations  also  suggest  that  Li7P3Sn  is  marginally  stable 
with  respect  to  decomposition  into  Li3PS4  and  Li4P2S7  in  the  sense 
that 

AH(Li7P3Sn)  =  AH(Li3PS4)  +  AH(Li4P2S7)  -  0.06  eV.  (5) 

Again,  this  is  within  the  expected  calculational  error. 

5.  Partial  densities  of  states 

The  partial  densities  of  states  can  be  evaluated  using  the  form 

Na(E)  s  e-(E-E"k>2^2.  (6) 

nk 

Here  the  smearing  factor  was  chosen  to  be  A  =  0.1  eV.  Wk  denotes 
the  Brillioun  zone  sampling  factor.  f“k  denotes  the  weighting  factor 
which  for  each  state  of  band  index  n  and  wave  vector  k,  was  taken 
to  be  the  corresponding  charge  within  the  atomic  sphere  at  site  a. 
The  sphere  radii  were  chosen  to  be  1.4bohr  for  O  and  S,  1.5bohr 
for  P  and  Cl,  and  1.6bohr  for  Li.  In  practice,  the  partial  densities  of 
state  for  each  atomic  site  are  averaged  over  all  similar  sites,  ( Na(E ) ). 
These  averaged  partial  densities  of  states  contain  qualitative  infor¬ 
mation  about  the  materials  in  terms  of  charge  transfers,  crystal  field 
splittings,  and  bond  formation.  In  all  of  the  phosphates,  the  bands 
are  dominated  by  the  2 p  states  of  O  while  for  the  thiophosphates, 
the  bands  are  dominated  by  the  3 p  bands  of  S. 

We  first  consider  materials  composed  of  isolated  phosphate  and 
thiophosphate  groups.  In  Fig.  5  the  partial  densities  of  states  of 
Li3PS4  and  Li3P04  are  compared.  The  S  3 p  valence  bands  of  Li3PS4 
cover  a  wider  range  of  energies  than  the  O  2 p  valence  bands  of 
Li3P04.  As  noted  in  our  earlier  work  [6]  the  two  lower  bands  have 
per  character  with  significant  hybridization  with  the  P  3s  and  3 p 
states,  while  the  top  of  the  valence  band  has  pn  character.  It  is 


E(eV) 


Fig.  6.  Partial  densities  of  states  for  Li7PS6  and  Li6PS5Cl.  Contributions  from 
unbounded  S  sites  are  indicated  with  “(u)”. 
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Fig.  7.  Partial  densities  of  states  for  U4P2S6  and  Li4P206. 


Fig.  8.  Partial  densities  of  states  for  U4P2S7  and  Li4P207.  Contributions  from  bridging 
S  or  0  sites  are  indicated  with  “(b)”. 

interesting  to  note  that  the  magnitudes  of  the  O  contributions  to 
(Na(E) )  the  partial  densities  of  states  are  generally  larger  than  those 
of  S,  indicating  that  the  charge  of  the  O  2 p  states  is  more  confined 
within  the  analysis  sphere  than  the  charge  of  the  S  3 p  states.  While 
the  calculated  band  gaps  within  density  functional  theory  are  sys¬ 
tematically  underestimated,  it  is  clear  that  the  gap  in  Li3PS4  is  less 
than  half  of  that  in  U3PO4. 


Fig.  10.  Energy  path  diagram  for  Li  ion  vacancy  migration  in  Li4P2S7  and  Li4P207.  The 
integer  labels  on  the  horizontal  axis  correspond  to  the  6  vacancy  sites  indicated  in  the 
structural  diagram  (Fig.  9).  Between  each  pair  of  optimized  vacancy  configurations 
are  7  NEB  image  configurations  along  the  minimum  energy  path. 

Next  we  consider  the  partial  densities  of  states  of  the  argy- 
rodite  structured  thiophosphates,  comparing  the  partial  densities 
of  states  of  Li7PS6  and  Li6PS5Cl  in  Fig.  6.  For  Li7PS6,  the  partial  den¬ 
sity  of  states  contribution  from  isolated  S  ions  (indicated  separately 
from  the  tetrahedrally  bonded  S’s)  are  found  to  contribute  rela¬ 
tively  narrow  bands  near  the  top  of  the  valence  band.  For  LigPSsCl, 
the  isolated  Cl  ion  contributions  are  at  lower  energy  than  those  of 
the  isolated  S  ions. 

In  Figs.  7  and  8  densities  of  states  for  Li4P206  and  Li4P2S6 
and  Li4P207  and  Li4P2S7  are  compared.  For  Li4P206  and  Li4P2S6, 
one  additional  valence  band  is  present  due  to  the  filled  P  3scr 
bond,  whose  contribution  is  concentrated  near  the  bottom  of  the 
spectrum.  In  contrast,  for  Li4P207  and  Li4P2S7,  while  states  with 
contributions  from  the  P  site  contribute  throughout  the  valence 
band,  there  is  no  additional  valence  state  due  to  filled  P  states.  On 
the  other  hand  there  are  two  types  of  0  (S)  states.  In  addition  to 
the  normal  tetrahedral  O  (S)  site  contributions,  the  “bridging”  O 
(S)  site  contributions  are  shown  separately  on  the  plot.  In  particu¬ 
lar,  there  is  a  band  at  the  bottom  of  the  valence  band  states  which 
corresponds  to  states  characterized  primarily  by  O  2 pa  or  S  3 per 
contributions  for  Li4P207  and  Li4P2S7,  respectively. 


Fig.  9.  Li  ion  vacancy  sites  used  in  migration  energy  study  in  Li4P2S7.  Similar  sites 
were  used  for  Li4P207. 


6.  Li  ion  migration  energies 

An  important  question  about  these  materials  is  how  the  Li  ion 
migration  differs  in  the  phosphate  and  thiophosphate  materials. 
The  ionic  conductivity  a  as  a  function  temperature  T  is  expected  to 
follow  the  Arrehenius  relation: 

a  T  =  Kne~Em/kT,  (7) 

where  k  is  the  Boltzmann  constant,  I<  is  a  material-dependent 
parameter,  n  is  the  concentration  of  mobile  Li  ions  or  Li  ion  vacan¬ 
cies,  and  Em  is  the  migration  energy. 

In  this  initial  study,  we  focus  on  the  migration  energy  Em  for  Li 
ion  vacancy  migration  in  idealized  crystals  of  Li4P207  and  Li4P2S7. 
As  discussed  in  Section  3,  these  two  materials  form  similar  stable 
and  meta-stable  crystals  with  P\  symmetry.  For  the  simulations, 
we  constructed  1x2x2  supercells  and  studied  the  migration  of 
a  Li  ion  vacancy  between  6  sites  separated  by  roughly  3.0  A  and 
3.6  A  for  Li4P207  and  Li4P2S7,  respectively.  The  sites  do  not  span  a 
complete  path  within  the  supercell  but  do  cover  a  representative 
range  of  crystal  environments.  The  6  sites  are  illustrated  in  Fig.  9 
for  the  Li4P2S7  supercell. 

The  band  occupancies  of  the  supercell  states  were  adjusted  in 
order  to  approximate  an  insulating  environment  for  ion  migration, 
and  the  excess  charge  was  compensated  by  adding  a  uniform  charge 
of  the  opposite  sign.  The  6  optimized  configurations  for  a  Li  ion 
vacancy  were  used  as  the  end  points  for  the  NEB  calculation  in 
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which  7  intermediate  images  were  used  to  approximate  the  min¬ 
imum  energy  path.  As  is  evident  from  the  energy  path  diagram 
shown  in  Fig.  10,  the  energy  landscape  for  this  system  is  com¬ 
plicated.  However,  the  results  indicate  a  clear  qualitative  result 
that  the  migration  energy  barriers  for  Li  ion  vacancy  migration  in 
Li4P2S7  is  less  than  half  that  in  Li4P207.  For  the  particular  migra¬ 
tion  path  chosen,  we  find  the  migration  energy  energy  barriers  to 
be  Em  =  0.3  eV  for  Li4P2S7  and  Em  =  0.8  eV  for  Li4P207. 

7.  Summary  and  conclusions 

This  study  has  found  some  interesting  similarities  and  differ¬ 
ences  between  the  phosphate  and  thiophosphate  materials.  Our 
structural  optimizations  were  generally  in  good  agreement  with 
the  experimental  structures,  although  there  were  some  obvious  dif¬ 
ferences.  For  example,  our  calculations  for  Li3PS4  indicated  that  the 
most  stable  structure  is  that  of  p-Li3P04,  while  the  experiment  [34] 
indicated  a  form  with  the  same  space  group  as  7-Li3  P04  with  partial 
occupancies  of  some  of  the  sites.  A  similar  discrepancy  was  found 
for  the  structure  of  Li4P2S6  [39].  These  differences  are  undoubtedly 
related  to  real  temperature  effects  which  are  not  considered  in  the 
simulations. 

In  general,  a  greater  variety  in  the  bonding  configurations  of 
the  thiophosphates  is  observed  compared  with  those  of  the  phos¬ 
phates.  The  argyrodite  structures  [35,36]  with  unbounded  S  and  Cl 
groups  within  the  crystal  and  increased  concentrations  of  Li  ions 
are  very  intriguing.  The  density  of  states  plots  shown  in  Fig.  6  indi¬ 
cate  that  the  unbounded  S  states  contribute  to  the  top  of  the  valence 
band.  The  corresponding  states  associated  with  unbounded  Cl  sites 
have  lower  energy  than  those  of  the  unbound  S  sites.  The  heats  of 
formation  for  these  materials  suggest  that  they  may  be  unstable  rel¬ 
ative  to  decomposition,  although  the  uncertainty  of  our  estimates 
do  not  allow  a  definitive  analysis. 

In  addition  to  our  analysis  of  the  stability  of  structures  and  sta¬ 
bilities,  we  have  presented  some  preliminary  results  pertaining  to 
Li  ion  vacancy  migration  within  a  representative  set  of  materials. 
We  find  that  the  energy  barriers  for  Li  ion  vacancy  migration  in 
Li4P2S7  to  be  less  than  half  that  for  Li4P207.  This  result  is  consistent 
with  the  experimental  observation  of  increased  ionic  conductivity 
in  the  thiophosphate  materials  compared  with  those  of  the  phos¬ 
phates.  Our  analysis  suggests  that  one  factor  contributing  to  the 
lowered  potential  barriers  in  the  thiophosphates  is  the  more  dif¬ 
fuse  valence  charge  distribution  near  the  S  sites  compared  with  the 
compact  valence  charge  distribution  near  the  O  sites. 

We  are  looking  forward  to  further  simulations  related  to  this 
intriguing  class  of  materials,  hopefully  with  collaboration  with 
experimental  studies. 
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